Quantum information processing with large nuclear spins in GaAs semiconductors 
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We propose an implementation for quantum information processing based on coherent manip- 
ulations of nuclear spins / = 3/2 in GaAs semiconductors. We describe theoretically an NMR 
method which involves multiphoton transitions and which exploits the non-equidistance of nuclear 
spin levels due to quadrupolar splittings. Starting from known spin anisotropics we derive effective 
Hamiltonians in a generalized rotating frame, valid for arbitrary /, which allow us to describe the 
non-perturbative time evolution of spin states generated by magnetic rf fields. We identify an ex- 
perimentally accessible regime where multiphoton Rabi oscillations are observable. In the nonlinear 
regime, we find Berry phase interference effects. 



PACS numbers: 76.60.-k, 42.65.-k, 03.67.-a 

Recent advances in spintronics (l) have shown that the 
coherent control of electron and nuclear spins in semi- 
conductors is experimentally feasible, enabling in partic- 
ular an all-optical NMR in GaAs, based on the hyperfine 
interaction between electrons and nuclei Such a 

control of nuclear spins can also be achieved via electri- 
cal gates as recently demonstrated for GaAs heterostruc- 
tures in the quantum Hall regime Q , or even via conven- 
tional NMR techniques directly accessing the nuclei ||. 
In the present work, we will show that such advances in 
coherent spin control have opened up the possibility to 
manipulate the nuclear spins / for the purpose of quan- 
tum information processing, thereby presenting a scheme 
that is based on ensembles of large spins I > 1/2 instead 
of qubits. Nuclear spins are ideal candidates for this pur- 
pose due to their long decoherence times. 

An implementation of the Grover algorithm || has re- 
cently been proposed for molecular magnets , based on 
a perturbative approach to the unitary Grover operations 
which encode and decode the information stored in the 
phases of small amplitudes a m ||]. An alternative ver- 
sion of Grover's algorithm was presented in Refs. |i||T(j| 
that is described by a Hamiltonian that lets a completely 
delocalized state of the form \ip) = Ylm=-i a m\ m ) > m 
some basis states |m) with equal occupation probabilities 
|am| 2 , evolve into a wanted localized state \M), where 
and \M) are degenerate and have a finite overlap. The 
information is encoded in the energies of |m). In order to 
produce \ip), we propose here a novel NMR scheme that 
allows us to generate any desired distribution of ampli- 
tudes a m , being not restricted to small values. For this we 
specifically exploit the properties of GaAs nuclei where 
quadrupolar spin splitting results in spin anisotropies and 
thus in non-equidistant energy levels- being a necessary 
condition for our scheme. The theoretical problem then 
is to find one magnetic rf pulse — inducing a unitary time 
evolution of the spins — that produces the desired spin 



state \yb) and a second rf pulse that lets \ip) evolve into 
\M), given certain spin anisotropies and adjustable mag- 
netic fields (see below). In a non-perturbative approach, 
we find an analytic solution to this problem, valid for ar- 
bitrary spin /. For the special case of GaAs with / = 3/2, 
we have confirmed our analytical results by exact numer- 
ics. In contrast to previous work Q| our method also 
holds for vanishing detuning energies, which turns out 
to be essential to perform non-perturbative unitary op- 
erations, i.e. quantum computations (QCs). Once the 
control over 21 magnetic fields is established, the scheme 
proposed here allows for QC and quantum storage with 
a single pulse, provided that there is sufficient signal am- 
plification due to the spin ensemble pd|. 
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FIG. 1. Multiphoton transition schemes for the coherent 
population of the I z eigenstates \m) of a nuclear spin / = 3/2. 
(a) Quantum computation (QC) scheme: The frequencies uik 
of the fields Hk are red (- •) and blue (- -) detuned. Diagrams 
containing blue detunings are negligible for large quadrupolar 
splitting, i.e. A ^> H8k > 0. (b) Rabi oscillation (RO) scheme: 
The magnetic fields H' k cos(tj' k t + 3>fc), k = 1, 2, 3, give rise to 
fc-photon RO. 
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As a first step towards this goal, it will be useful to gen- 
erate and monitor multiphoton Rabi oscillations, as we 
describe in detail below. Finally, we show that oscillating 
quadratic transverse spin terms, which can be generated 
by optical pulses in GaAs PJl^] , give rise to Berry phase 
oscillations [H in the transition probabilities. 

In the following we mainly focus on a nuclear spin of 
length I = 3/2, as appropriate for GaAs, but indicate its 
generalization to arbitray I. Our nuclear spin system is 
described by the Hamiltonian Ho = Hz + 7Yq, consist- 
ing of the nuclear Zeeman energy Hz = —gNfJ-NH z I z , 
gN = 1-3 p|, and the quadrupolar splitting |5|] Hq — 
A[3I Z — 1(1 + 1)]. The quadrupolar constant is A = 
7 x KT 7 K for 69 Ga, A = 3 x KT 7 K for 71 Ga, and 
A = 2 x KT 6 K for 75 As nuclei §. For the purpose 
of QC we need to achieve complete control over unitary 
state evolutions, i.e. control over amplitudes a m to form 
a desired superposition — Yl m =-I 0,7,1 l m ) °^ ^ ne nu ~ 
clear basis states |m) (eigenstates of Ho)- Our goal is 
now to show that such a control over a m is indeed feasi- 
ble under experimentally attainable conditions. 




FIG. 2. Preparation of js) = (l/V%) l m ) h y 

means of Eq. (2) in the QC scheme, which takes about 0.2 
ms for Hi = H 2 = 1 G, H 3 = 0, Si = 6083 s"\ and 5 2 = 0. 
The duration of the QC operation is < l/2^ a 'bi- The analyt- 
ical result is confirmed by numerics. 

We start from a configuration where mainly the ground 
state |3/2) is populated, see Fig. |Tj. This can be achieved 
by the Overhauser effect The next goal is to co- 

herently populate all or a part of the excited states 
|m), to 7^ 3/2, by means of Am = 1,2,3-photon tran- 
sitions. Fig. [I] shows the two transition schemes, QC 
and RO, which will turn out to be appropriate for quan- 
tum computation (QC) and multiphoton Rabi oscilla- 
tions (RO), resp. In the QC scheme the frequencies 
LUk of the external transverse magnetic fields, H X:k (t) = 
H k (t) cos(u k t + k = 1,2,3, are blue (5 k < 0) and 
red (5 k > 0) detuned with respect to the energy level 
separations tiu! mm >. In the RO scheme, the transverse 
fields H' x k (t) = H' k (t)cos(u' k t+® k ), k = 1,2,3, oscil- 
late at frequencies uj' Am = ^3/2-Am.3/2/Am, which are 
blue detuned by 3 A (6 A) for the two(three)-photon tran- 



sition. For GaAs, Wfc,u^ ~ 10 MHz with S k ~ 1 kHz, 
and a longitudinal magnetic field H z ~ 1 T is appro- 
priate. It is desirable to make H z sufficiently large to 
accommodate many spin precessions before the spins de- 
phase. We note that in contrast to the fields H Xj k(t), 
the fields H' x k (t) lead to transitions governed by non- 
commuting operators, with the important consequence 
that the RO scheme suffers from strong interferences 
between the transitions if two or more fields H' x k (t) 
are nonzero, leading to a quick loss of amplitude con- 
trol. Indeed, the RO scheme allows control of a m 's only 
for times t < 2K(V k ' + V k ,)/V k V k ,, which we estimate 
from the Baker-Campbell-Hausdorff formula and which 
we confirmed by exact numerical calculations. Here, 



2[(gNVNH k ) k p 3/2 _ k ^ /2 ]/l\ k - = ±hLJz_ :j ^ (see be- 
low). Although the RO scheme is only suited for QCs 
using perturbative approaches, it has its usefulness for 
testing the coherence of the spin system (see below) . 
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FIG. 3. Grover algorithm calculated by means of 
Eq. (2) in the QC scheme (numerically confirmed), where 
\s) — (l/y/3) X/m=-i/2 l m ) ^ s concentrated mainly into 
| — 1/2} after 0.55 ms for H2 = h52/2gN/-iN = 1 G, hi = h<2, 
h 3 = 0, Si = 0. The duration of the QC is < l/2f| 



/2) 
'Rabi 



Now we proceed with demonstrating the existence of 
the desired spin transitions in the QC scheme. For this 
we evaluate the transition amplitudes for the diagrams 
of Fig. 1 (a) in high-order perturbation theory which al- 
lows us then to obtain an appropriate non-perturbative 
Hamiltonian (see below). The three transverse fields 
H x ,k(t) complete the Hamiltonian H = Ho + V(t), 
where V(t) = £)a=i 5jvMiV-fffc(*) cos(cJfc£ + $ k )I x , with 
I x = (J + + I-)/2, and phases & k (see below). Then 
we expand the S'-operator S = SJlo ^ m P° wers of 
V(t). We use rectangular pulse shapes of duration T for 
all fields, i.e. H k (i) = H k for -T/2 < t < T/2, and 
otherwise. Then we obtain 
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Efe=l W fe)i Pm,m' = nr=m( /s l / -l fc + 1 >' and <^ (t) m = 
T- f + X 7 9 2 e iuJt dt = sin( " T/2) is the delta-function of width 

Z7T J —1 / 2 7TUJ 

1/T. The energy is conserved for wT ^> 1. Also, 
the duration T of the rf pulses must not exceed the 
dephasing time of the spin states. Interestingly, 

(3) (2) 
lim,4_ ) .o >S_3 / 2 3 / 2 = hmA-fO <S_i/ 2 3/2 = 0> i- e - destruc- 
tive interference is maximal. However, if A 3> 7i|<5fc|, 
k = 1,2,3, destructive interference is negligible. 
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FIG. 4. Numerical solution for the three-photon ROs of 
75 As nuclei between |3/2) and |— 3/2) driven by H' A = 20 G 
with the RO scheme (b) in cw mode. H[ — H' 2 = 0. 

Now we are in the position to extract an effective 
Hamiltonian, which governs the desired unitary evolu- 
tions §,[l(J. For the QC scheme we use the Hamiltonian 
H.. After applying the rotating wave approximation |^] 
we keep only the most left diagram of Fig. 1(a), which 
gives the dominant contribution to the transition ampli- 
tudes for h\Sk\ <C \A\. This is a direct consequence of 
the non-equidistance of the energy levels \m) due to the 
quadrupolar splitting. It is now possible to eliminate the 
time-dependence of H. by a unitary operation U(t), the 
matrix elements of which can be determined by solving 
21 linear equations. This is a transformation to a gen- 
eralized rotating frame. Then, for a spin / we obtain an 
effective time-independent Hamiltonian 
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where h k = g N fj, N H ky /k(2I + I - k)/2 (k = 1,...,2I). 
Focusing on / = 3/2, we obtain e.g. for H3 — ap- 
proximately = Ho + 
hie i{ Ul t+^) | 3 /2) (l/2| + /j 2e i(^t+* 2 ) |i/ 2 ) (_l/2| +h.c. 
Applying U(t) = e -<[(«i+« 2 )t+(»i+*B)]/2 13/ 2 ) (3/2| + 

e t[(u;i-W2)t+(*i-*a)]/2 1 2/2) (1/2| + 

el [( Wl +^)t+(*i+*2)]/2|„ 1 /2)(-i/2| yields «PJ. Note 
that the Hamiltonian in Eq. (^) remains valid even in 
the limit 5k —> 0, where perturbation expansions such as 
in Eq. (|l|) break down. However, we must require that 
\gN^NHk\ <C I .A I, which means that the larger \A\, the 

faster the QCs. Propagators of the form W (t)e~ lH '<> t>t / h 
have 21 phases $k and 21 detunings H5k, which deter- 
mine the 21 phases and the 21 moduli of a m |16|. 

For Grover's algorithm [p|JlO|] we must first produce 
\s) = (l/\/n)E ra l m ) ( see Fig. H): n being the number 
of basis states involved in the search. Then we make 
use of the degeneracy between \s) and |M), which yields 
the resonance condition h k = %§i—m/'^ / Vfc, if 5k = 
Vfc ^ I - M. In contrast to ft^i? has only 

nearest-neighbor coupling, which results in a decreasing 
amplification of \M) with increasing / or \M\. However, 
even for the largest nuclear spin / = 9/2, we find that the 
resolution for identifying \M) is still sufficient (> 10%). 
Fig. H shows the example where \M = —1/2) is found out 
of the three states \m), m = 3/2, 1/2, -1/2, for I = 3/2. 




FIG. 5. Numerical solution for the two-photon ROs of As 
nuclei between 1 3/2) and | — 1/2), driven by H' 2 = 10 G accord- 
ing to the RO scheme in cw mode, and H[ = H3 — 0. 

As a first test for the proposed schemes, it would be 

useful to measure generalized ROs involving multiphoton 

absorptions. They can be thought of as nutation of the 

large spin / between spin states \m). First, we consider 

the QC scheme. For the two-photon RO, with frequency 
(2) 

^Rabi' t° become observable, we need \gN^NHk\ <C Mi, 

k = 1,2, so that the one-photon transitions arc com- 

(2) 

pletely suppressed. To obtain v^ hi , it is useful to think 
of (0) as describing the dynamics of a (fictitious) particle 
in a triple well with nearest-neighbor tunnel coupling hk ■ 
The independent control of the tunnel couplings hk and 
the biases h5k between the wells is ensured by a large 
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value of A. Then, the energy ("tunnel") splitting fTs] ] 
between |3/2) and | — 1/2) reads for 82 — 



i (2) 

i Rabi 



V3(g Nf x N ) 2 H 1 H 2 /S 1 



(3) 



("2") (2) 

which gives ^p abi = Aj^ bi /27r?i. In order to obtain large 

(2) 

Rabi frequencies fj[ abi , the external fields H\,H2 and 
the detuning fi8\ are to be maximized under the condi- 
tions \Hx\, \H 2 \ < ^l^l/ffjvMiv < |A|/ffArAtAr @, i.e. the 

(2) — 
larger |A| the larger fjtabi can be achieved. We note that 

|A| could e.g. be enhanced by optical laser pumping ||] 
or by modulated electric field gradients p"^ ]. 

Next we turn to the RO scheme. Here it is sufficient to 
apply only one single field H' x k (t) in order to see the mul- 
tiphoton ROs shown in Figs, jlj, || We now also allow for 
oscillating quadratic transverse anisotropics which can 
be externally generated by modulating the electric field 



gradient felt by the nuclei 
most general Hamiltonian [5| 



For this we adopt the 



H' = A[3I 2 - 1(1 + 1)] - g Ni i N H z I z + e l < th [-h' k I x 



B{I x I z +I z I x )+C(l 2 -I 2 



-iuj' k tl z 



(4) 



where h' k = g^^NH k (k = 1,2 or 3). Next we trans- 
form TC to the rotating frame, which yields 7i' ro t = 
A[3I 2 -I(I+l)}-(g N ^ N H z I z -hcu' k )I z +B (I X I Z + I Z I X ) + 
C (I 2 — I 2 ) + h' k I x . Then the time evolution takes the 
simple form \ij)(t)) = e ^' k ti z e -iH' rat t/h ^ Although the 
transverse quadratic term C is not in resonance with any 
transition energy, it leads to a time-independent trans- 
verse quadratic anisotropy in the rotating frame. For 
the 3-photon transition in the RO scheme we obtain the 
following Hamiltonian in the rotating frame 



o//(3) _ 
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where we have neglected the B term since we choose 

B < h' 3 . Inserting a typical value C = -10 -10 K |,we 

(3} 

obtain oscillations of the splitting Agi between |3/2) 
and |— 3/2) as a function of H^, see Fig. ^. These oscilla- 
tions are due to the Berry phase in a biaxial spin system 
as shown in Note that C must be negative for the 
Berry phase interference to occur [O. Also, the Berry 
phase interference is present only for Am-photon transi- 
tions with Am > 2. In Figs. [|and [|the population prob- 
abilities |a m (t)| 2 are shown for C = 0. The corresponding 
normalized magnetization reads M(t) = ^2 m m\a m \ . 

In order to build parallel single-spin quantum comput- 
ers or high-density memory devices, one could apply an 
inhomogeneous magnetic field to the GaAs sample. Then 
the nuclear Larmor frequencies are spread over a wide 
range, which could be divided into small frequency inter- 
vals that are individually accessible by NMR sources. 



In conclusion we have shown that via multiphoton 
transitions a controlled superposition of spin states can 
be achieved by appropriate field pulses. Our results can 
be extended to arbitrary spin / and to any single-particle 
quantum system with non-equidistant energy levels. 
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FIG. 6. Berry phase oscillation. The three-photon transi- 
tion probability vanishes where A'^2 bi is zero. 
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